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Abstract. Gamma-ray emission from classical novae is 
dominated, during the first hours, by positron annihilation 
resulting from the beta decay of radioactive nuclei. The 
main contribution comes from the decay of 18 F and hence 
is directly related to 18 F formation during the outburst. A 
good knowledge of the nuclear reaction rates of production 
and destruction of 18 F is required to study 18 F synthesis 
in novae and the resulting gamma-ray emission. The rates 
relevant for the main mode of 18 F destruction (i.e, through 
proton captures) have been the object of many recent ex- 
periments. However, subsequent analyses were focused on 
providing rates for X-ray burst nucleosynthesis not valid 
at nova temperatures (lower than 3.5xl0 8 K). Accord- 
ingly, it is crucial to propose and discuss new reaction 
rates, incorporating all new experimental results, down to 
the domain of nova nucleosynthesis. We show that in this 
temperature regime, the 18 F(p,7) 19 Ne and 18 F(p,a) 15 
reaction rates remain uncertain and deserve further exper- 
imental and theoretical efforts. Our hydrodynamic calcu- 
lations including the new nuclear rates demonstrate that 
their impact on 18 F synthesis in nova explosions is quite 
large and, consequently, the early gamma-ray emission 
from classical novae is also affected. 

Key words: Nuclear reactions, nucleosynthesis, abun- 
dances - Novae, cataclysmic variables - Gamma rays: the- 
ory 



1. Introduction 

Classical novae emit gamma-ray radiation at and below 
511 keV during the early epochs after the explosion. This 
emission is produced by electron-positron annihilation in 
the expanding envelope, and the subsequent Comptoniza- 
tion of the resulting gamma-ray photons, and it shows a 
line, at 511 keV, and a continuum, between 20 and 511 keV 
( pomez-Gomar et al. 1998 ) . The positrons responsible for 
t his emission come main l y from the disintegration o f 18 F 
( Leising fc Clayton 1987 ; Gomez-Gomar et al. 1998 ), be- 
cause its lifetime (r— 158 min) is such that positrons are 



emitted at the "right time", i.e., when the expanding en- 
velope starts to be transparent to gamma-ray radiation. 
Therefore, the amount of radiation emitted strongly de- 
pends on the 18 F content of the nova envelope. 

The synthesis of 18 F in novae depends largely on 
some key nuclear reaction rates of 18 F destruction 
and production which are far from being well known. 
This is the case, in particular, of the 18 F(p,7) 19 Ne 
and 18 F(p,a) 15 reactions. Recent experimental stud- 
ies ( praulich et al. 1997 ; Utku et al. 1998| ) drastically im- 
proved the knowledge of these reaction rates with respect 
to previous studies ( Wiescher fc Kettncr 1982| ). In a re- 
cent paper ( Hernanz et al. 199S), we ha ve analyzed the 
influence of these rates ([Utku et al. 1998|) of 18 F(p,7) 19 Ne 
and 18 F(p,a) 15 on the final yields of 18 F for different 
models of CO and ONe novae. The effect of the new rates 
was important, since a factor of 10 reduction in the yields 
and in the resulting gamma-ray fluxes was obtained for all 
the models. Therefore, we concluded that a more detailed 
analysis of the reaction rates was necessary, in order to 
predict the gamma-ray emission of classical novae. 

The rates proposed by Utku et al. (1998) were limited 
to relatively high resonance energies and temperature do- 
mains, more appropriate for temperatures typical of X- 
ray bursts than for those of classical novae. The reason is 
that for their determination neither the influence of un- 
certainties on low energy resonance strengths (based on 
assumed reduced width) nor the effect of all low energy 
resonance tails were considered. At higher temperatures 
the rates are more reliable since their main contributions 
come from two directly measured resonances. The pur- 
pose of this paper is to provide rates for the 18 F(p,7) 19 Ne 
and 18 F(p,a) 15 reactions valid in the domain of tem- 
perature of nova nucleosynthesis, incorporating the latest 
experimental data. In addition to the nominal rate, we 
provide upper and lower limits. Since other nuclear reac- 
tions also affect 18 F synthesis in novae, a global analysis 
is done , including results f rom the recent NACRE compi- 
lation ( Angulo et al. 1999 ). With the new rates affecting 
18 F synthesis, new nova models have been computed, in 
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order to determine the mass of 18 F they eject and the im- 
pact of the new yields on their early gamma-ray emission. 

The organisation of this paper is the following. In sec- 
tion H we describe the nucleosynthesis of 18 F. In section || 
we discuss in detail the 18 F+p rates and the correspond- 
ing uncertainties while in the following section (||), we 
briefly discuss other recently published rates. In section 
H we present new results that show the influence of these 
new rates on 18 F production. Following the conclusion, 
an appendix gives some analytical approximations to the 
rates. 



2. Synthesis of 18 F in classical novae 

In this Section, we will focus on the nuclear paths leading 
to 18 F synthesis, through a detailed analysis of a repre- 
sentative model of an ONe nova : a 1.25 M Q white dwarf, 
accreting solar-like matter at a rate of 2 x 1CP 10 M Q .yr _1 
, assuming a 50% degree of mixing with the outermost 
core of composition taken from Ritossa, Garcfa-Berro & 
Iben (1996) (mainly 16 and 20 Ne with traces of other Ne, 
Na and Mg isotopes). The evolutionary sequences leading 
to a nova outburst have been computed with an updated 
version of the code SHIVA (see Jose & Hernanz 1998), 
a one-dimensional, implicit, hydrodynamical code in La- 
grangian formulation, which follows the course of the ex- 
plosion from the onset of accretion up to the expansion 
and ejection stages. The nuclear reaction network is de- 
scribed in Jose et al. (1999), and in Hernanz et al. (1999). 
In particular, we use the reaction rates for proton captures 
on 18 F and 17 based on Utku et al. (1998) and Landre 
et al. (1989), respectively. 

The synthesis of 18 F in nova outbursts takes place 
within the hot-CNO cycle (see fig. [l]). Regardless of the 
nova type (either CO or ONe, according to the com- 
position of the underlying white dwarf), the initial 16 
abundance is quite large (up to 25% by mass). Hence, 
initial 16 is the main source of the formation of 18 F, 
either through the chain 16 0(p,7) 17 F(p,7) 18 Ne(/?+) 18 F 
or through 16 0(p,7) 17 F(/3+) 17 0(p,7) 18 F, which reflects a 
competition between (p,7) reactions and /3 + -decays. 

Snapshots of the evolution of isotopes relevant to 18 F 
synthesis (i.e., 16 < 17 0, 17 ' 18 > 19 F and 18 > 19 Ne) are shown 
in Fig. At the early stages of the explosion, when 
the temperature achieved at the burning shell reaches 
Tb s ~ 5 x 10 7 K (Fig. [2], first panel), the main nuclear 
activity in this mass region is driven by 16 0(p,7) 17 F, fol- 
lowed by its /3 + -decay to 17 0. However, the temperature 
is not high enough to burn a significant fraction of 16 O. 
On the contrary, the abundance of 17 rises by several 
orders of magnitude with respect to its initial content. 
Both 17 ' 18 F are being synthesized at the burning shell by 
means of proton captures on 16 ' 17 0, respectively, whereas 
the amount of 19 F remains essentially unchanged. The mi- 
nor activity driven by proton capture reactions onto 17 > 18 F 



results also in a moderate increase of both 18,19 Ne (below 
1CP 10 by mass). 

A similar trend is found when Tj, s ~ 7 x 10 7 K (Fig. Q 
second panel). The temperature rise increases the num- 
ber of proton captures onto 16 0, leading to 17 F, which 
rapidly decays into 17 0. The effect of convection, which 
already extends through most of the envelope (i.e., ~ 150 
km), is shown in the smooth distribution of both 17 > 18 F, 
previously synthesized in the burning shell and efficiently 
carried away to the outer envelope. 19 F is in turn re- 
duced down to 50% due to 19 F(p,a) 16 0, which dominates 
19 Ne(/3+) 19 F. 

When T bs reaches 10 8 K (Fig. [|, third panel), 17 
shows a flat profile along the envelope, with a mean 
mass fraction of - 3 x 10~ 3 , due to 16 0(p,7) 17 F(/3+) 17 0, 
which dominates destruction through 17 0(p,7) 18 F. With 
respect to fluorine isotopes, we stress an important syn- 
thesis of 17 ' 18 F at this stage, driven by 16 0(p,7) 17 F 
and 17 0(p,7) 18 F, whereas 19 F is almost fully destroyed. 
18,19 Ne continue to rise due to proton captures onto 17 ' 18 F. 

A dramatic change is found as soon as the burning shell 
reaches 2 x 10 8 K (Fig. @, fourth panel). At this time, the 
mass fraction of 16 has been reduced down to ~ 3 x 10~ 2 . 
For the first time, 17 decreases (by one order of mag- 
nitude with respect to the amount shown in previous 
panel), since destruction through 17 0(p,7) 18 F dominates 
17 F(/3 + ) 17 in the vicinity of the burning shell. On the 
contrary, 17 F exhibits a significant rise, up to ~ 3 x 10 -2 
by mass (i.e., one order of magnitude higher than in previ- 
ous panel), so that the amount of 17 F becomes larger than 
that of 17 in almost all the envelope. Again, 16 0(p,7) 17 F 
dominates both 17 0(p,7) 18 Ne and 17 F(/3+) 17 near the 
burning shell. The evolution of 19 F is still dominated by 
19 F(p,a) 16 0, which in fact has no significant influence on 
the 16 O content. Also noticeable is the dramatic rise of 
both 18,19 Ne, which increase by several orders of magni- 
tude, mainly due to the fact that at such temperatures, 
proton captures onto 17 > 18 F become faster than the cor- 
responding /3 + -decays. It is worth noticing that at this 
stage, a non-negligible path leading to 18 F synthesis is 
now driven by 18 Ne(/3 + ) 18 F, basically at the outer en- 
velope shells, which compensates an efficient destruction 
through proton captures near the burning shell. 

The burning shell achieves a maximum temperature of 
2.51 x 10 8 K (Fig. [| fifth panel). 16 is reduced to 0.096 
by mass in the burning shell. Whereas 17 is destroyed by 
(p,a) and (p,7) reactions at the burning shell, its mean 
abundance in the envelope increases due to 17 F(/3 + ) 17 
which dominates in almost all the envelope. On the con- 
trary, 17 F is generated in the burning shell by means of 
16 0(p,7) 17 F, whereas it is being destroyed by /3+-dccays in 
the outer envelope. The evolution of 18 F reflects a compe- 
tition between destruction and creation modes at different 
locations of the envelope: whereas its content is reduced 
by both 18 F(p,a) 15 reactions at the burning shell and by 
18 F(/3 + ) 18 at the outer envelope, a dominant source for 
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18 F synthesis through 17 0(p,7) 18 F is found at the inter- 
mediate shells. Proton captures onto 17 ' 18 F and convective 
transport continue to pump 18 ' 19 Ne to the outer envelope 
shells, which is at the origin of the rise of 19 F at this stage 
(through 19 Ne(/3+) 19 F). 

Shortly after, due to the sudden release of energy from 
the short-lived species 13 N, 14 > 15 and 17 F, the envelope 
begins to expand. As a result of the drop in temperature, 
(p,7) and (p,a) reactions are basically restricted to the 
vicinity of the burning shell, whereas most of the enve- 
lope is dominated by /3+-decays (Fig. sixth and seventh 
panels). Hence, whereas 17 is powered by 17 F(/3+) 17 
along the envelope, the amount of 17 F and 18 ' 19 Ne de- 
creases significantly. The 18 F abundance at these late-time 
stages of the outburst increases due to 17 0(p,7) 18 F, which 
is dominant at the intermediate layers of the envelope 
(whereas in the vicinity of the burning shell 18 F is effi- 
ciently destroyed by 18 F(p,a) 15 0, the outer envelope is 
dominated by 18 F(/3+) 18 0). Moreover, destruction of 19 F 
through 19 F(p,a) 16 at the bur ning shell is dominated by 
19 Ne(/3+) 19 F elsewhere. 

The final stages of the outburst (Fig. Q eighth panel), 
as the envelope expands and cools down, are dominated 
by the release of nuclear energy by /3+-decays such as 
18 F(/?+) 18 0, or 17 F(/3+) 17 0. The resulting mean abun- 
dance of 18 F in the ejected shells in this Model is X( 18 F) = 
2.7 x 10~ 4 . Most of the envelope is, however, dominated 
by the large abundances of 16 < 17 (X( 16 0)= 6.8 x 1CT 2 , 
X( 17 0)= 3.9 x 1CT 2 ). A residual 17 F, which is still decay- 
ing into 17 0, and a non-negligible amount of 19 F (3 x 10~ 6 , 
by mass) are also present. 

In summary, the synthesis of 18 F in classical novae 
is essentially controlled by five proton-capture reactions, 
16 0(p, 7 ) 17 F, 18 F(p, 7 ) 19 Ne, 18 F(p,a) 15 0, 17 0(p, 7 ) 18 F, 
17 0(p,a) 14 N, and several /3+-decays. Accordingly, the cor- 
responding reaction rates deserve further attention and 
their uncertainties are discussed in Sect. || for the 18 F+p 
rates and Sect. ^ for the other rates. 

3. The 18 F(p,a) 15 and 18 F(p,7) 19 Ne reaction 
rates 

Uncertainties on the 18 F(p,7) 19 Ne and 18 F(p,a) 15 reac- 
tion rates used to reach many orders of magnitude due to 
the very limited spectroscopic data available for 19 Ne in 
the domain of interest. The rate estimated by Wieschcr 
& Kettner (1982) has become obsolete by the recent mea- 
surements of Rehm et al. (1995), Coszach et al. (1995), 
Rehm et al. (1996, 1997) and in particular by the work of 
Graulich et al. (1997), Utku et al. (1998) and Butt et al. 
(1998). We concentrate on the faster 18 F(p,a) 15 reaction 
but most of the discussion also applies to the 18 F(p, 7 ) 19 Ne 
one. 

Before 1997 and below the excitation energy (E x ) cor- 
responding to the resonance energy (E r ) of 0.4 MeV, only 
two levels were known E x = 6.437 and 6.742 MeV but 



with unknown or uncertain spins and parity (J^)^} On 
the contrary, various levels were known in the correspond- 
ing region of 19 F, the conjugate nucleus (see Fig. [||). For 
instance, in the E r =0.—1. MeV region, only 6 levels were 
known and many spins unknown in 19 Ne, while 19 levels 
had been observed in the corresponding 19 F region. Hence, 
for the estimation of the rates, Wiescher & Kettner (1982) 
considered three known levels at E x = 6.437, 6.742 and 
6.862 MeV and postulated three others. The correspond- 
ing rate was highly uncertain since the location of these 
three postulated levels was approximate and because other 
levels found in 19 F are expected to have counterparts in 
19 Ne in the region of interest. This is the case, in partic- 
ular, of the E x , r = 6.429, 1/2", 6.497, 3/2+ and 6.527, 
3/2+ levels of 19 F corresponding to strong expected res- 
onances because of their low centrifugal barrier, i.e. low 
transfered orbital angular momentum (£ p ) of 1, and 
respectively. 

Even though the recent experiments have not been able 
to find all the counterparts of the 19 F levels (see Fig. [||), 
they provided direct measurement for the strengths of two 
resonances which likely dominate the rate in the domain 
of nova nucleosynthesis and the location of several new 
19 Ne levels. The strength of the resonance correspond- 
ing to the E x = 6.742 MeV level (£V=330 keV) has re- 
cently been measured directly by Graulich et al. (1997) 
with a 18 F beam provided by the Louvain-La-Neuve fa- 
cility. This level is thought to be the analog of the 6.787, 
3/2~ 19 F level. The measured strength (ojj), 3.5 ± 1.6 eV 
( praulich et al. 1997Q , is in good agreement with the Wi- 
escher & Kettner (1982) estimate and accordingly does not 
induce a significant change in the rate. On the contrary, 
the resonance associated with the E Xl J^ = 7.067 MeV, 
3/2 + level (-E r =659 keV) strongly modifies the rates. It 
has been studied by Rehm et al. (1995, 1996), Coszach 
et al. (1995), Graulich et al. (1997), Utku et al. (1998) 
and Butt et al. (1998) and its s trength has been me asured 
directly ( Coszach et al. 1995[ |Graulich et al. 1997 ). It is 
located well outside of the Gamow peak but due to its 
total width (r «30 keV) , the contribution of its tail alone 
is greater than the Wiescher & Kettner (1982) rate in the 
domain of nova nucleosynthesis. However, no counterpart 
of this broad level was known in 19 F. The closest 3/2+ 
19 F level is located at E x = 7.262 but its width is much 
smaller (r <6 keV). To remove this incompatibility, Rehm 
et al. (1996) claimed that the 7.067 MeV width could be 
smaller (r «14 keV). A 7.114, 7/2+ level was known in 
19 F with the possibility that it hides an unresolved J 71 = 
3/2 + level as suggested by Smotrich et al. (1961). A re- 
cent gamma ray spectroscopy study ( Butt et al. 1998| ) of 
19 F produced by the 15 N(cv,7) 19 F reaction has confirmed 
the presence of a 7.101 MeV, 3/2+, T = 28 ± 1 keV level. 



1 For more detailed explanations of the standard nuclear 
physics notations used, see e.g. the appendix in Jose et al. 
(1999). 



4 



Alain Coc et al.: Influence of new reaction rates on F production in novae 



Hence the high width found by Coszach et al. (1995) and 
Utku et al. (1998) for the E X ,F = 7.067 MeV, 3/2+, 19 Ne 
level can be understood. (Note however that this analog 
assignment has been questioned very recently by Fortune 
& Sherr (2000).) 

Utku et al. (1998) have also found three new 19 Ne lev- 
els and have made tentative assignments of analog lev- 
els in 19 F (see FigjJ). They recognized the 6.437 MeV, 
1/2" level (E r =2& keV) as the analog of the 6.429, 1/2", 
r = 280 keV one in 19 F from its large measured width 
(216 ± 19 keV). Even though it is located at a very low 
energy (in the context of novae) it is so broad that its 
tail may lead to a significant contribution to the rate (see 
Fig.[|). The 6.449 MeV 3/2+ (£ r =38 keV) level is one 
of the three new ones found by Utku et al. (1998). It is 
also of great importance since the contribution of its tail 
can dominate the astrophysical factor in the relevant en- 
ergy range (see Fig. [J). In comparison the two other new 
levels give a smaller (6.698 MeV 5/2+, £ r =287 keV) or 
negligible contribution (6.419 MeV 3/2+, E r =8 keV). 

A new rate has been provided by Utku et al. (1998), 
using their new data together with available spectroscopic 
information on 19 Ne and 19 F supplemented by estimates 
of radiative and proton widths when missing. This rate 
is already higher by a factor of up to « 3 x 10 2 when 
compared with Wiescher & Kettner (1982). Nevertheless, 
even though this rate is now set on a firmer basis than 
the Wiescher & Kettner (1982) one, it remains quite un- 
certain as it depends directly on the estimated proton 
widths. Except for the two resonances whose strengths 
have been measured directly by Graulich et al. (1997), 
the proton widths are assumed to be equal to 0.1 or 0.01 
Wi gner limit, i.e. 9^= 0.1 or 0.01, according to their par- 
ity ( Utku et al. 1998 ). For the resonances considered here, 
one has T p <C Txot ~ Tq so that their strengths are 
(ujj ~ ujT p ) proportional to chosen # 2 values. Accordingly, 
to obtain lower and upper limits for the rate, we used the 
two extreme cases 9^=0 and 9^=1 respectively and re- 
gardless of parity but kept Utku et al. (1998) # 2 values for 
the nominal rate (see Fig.[6|). The lower limit is quite cer- 
tain as it is given by the contribution of the two directly 
measured resonances whose parameters are now reliably 
determined. It remains higher than the Wiescher & Ket- 
tner (1982) rate above ss50xl0 6 K. The nominal rate is 
conditioned by the rather arbitrary choice of 9^=0.1 or 
0.01 that we adopt following Utku et al. (1998). However, 
our nominal rate is significantly higher [up to a factor 
of «10 and 30 around 10 8 K for (p,a) and (p,7) reac- 
tions] than the one from Utku et al. (1998) for typical 
nova temperatures (see Fig.^BJ) . This is due to the inclu- 
sion of the contribution of the tail (Ftm = r Q =4.3 keV) of 
the 6.449 MeV level (see Fig. [4|) not considered by Utku 
et al. (1998). The corresponding rates (p,a) are presented 
in FigJ^ and have been calculated by numerical integra- 
tion of the Breit- Wigner formula for all broad resonances. 
For the (p,7) rates, we used the same procedure to ob- 



tain the low, nominal and high rates (see Figs. [5] and Q). 
Radiative widths were taken from 19 F analog levels as in 
Utku et al. (1998) except for the important resonance at 
659 keV as Bardayan et al. (1999) provided an experi- 
mental value (r 7 =0.39 eV) for the analog level. We also 
adopted the direct capture (DC) contribution from Utku 
et al. (1998) even though the spectroscopic factors used 
( 18 0®p instead of 18 F<g)p or even 18 F<8>n) may not be the 
more appropriates. These rates have been calculated using 
all the available spectroscopic data collected in Table II of 
Utku et al. (1998) except for the inclusion of a new gamma 
width for the 659 keV resonance and for a different proton 
width for the 26 keV resonance [1% of the Wigner limits 
i.e. 2.5xl0 -17 eV, as for other negative parity resonances, 
instead of 6.6xl0~ 20 eV of unknown origin in Utku et 
al. (1998)]. These numbers are taken at face value to cal- 
culate the nominal rate only. For the low and high rates 
unknown proton widths are allowed to vary between and 
the Wigner limit as discussed above. It is clear that this 
range (0 < 9^ < 1) represents the most extreme values as 
are the corresponding rates. This should be kept in mind 
when interpreting the astrophysical implications. 

However a few points have been neglected due to the 
lack of experimental or theoretical information. First of 
all, the identification of all E x 6.6 MeV analog levels is 
not complete yet (see Fig. [^) and is not certain. In particu- 
lar, new contributions to the nominal and high rates corre- 
sponding to missing levels cannot be ruled out. The tail of 
the 38 keV resonance has been calculated using the large 
total width derived from the 19 F analog level proposed 
by Utku et al. (1998). Another 3/2+ level, with smaller 
width, lies 30 keV below in 19 F. It has been assigned 



(Utku et al. 1998) to be the analog level corresponding 
to the 8 keV (E.^6.419 MeV) resonance in 18 F+p (see 
Fig. [^) . It is possible that the analog assignment for these 
two levels is inversed or that they are mixed but at least 
one or the other resonance (8 or 38 keV) should be broad 
enough to dominate the rates at low nova temperature. 
For the (p,7) rates, we neglected any difference of radiative 
widths between analog levels and the uncertainty on the 
DC contribution in front of the uncertainty on the proton 
widths. These rates do not either include possible interfer- 
ence effects arising from the interference between the 3/2+ 
resonances located at 38 and 659 keV. In a favourable case 
(when 9^ «0.01 for the 38 keV resonance), destructive in- 
terference could reduce the low rate by a factor of «10 in 
the sensitive T « 10 s K region. However, at present the 
rates provided in the appendix are those we recommend. 



4. Other reactions affecting 18 F production 

According to the recent compilation of charged-particle 
induced thermonuclear reaction rates NACRE (Angulo 
et al. 1999), the rate for 16 0(p,7) 17 F is rather well 
known and suffers from a little uncertainty (i.e., a 
factor of ~ 2 between the high and low rates). Even 
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though 17 F has a much shorter lifetime than 18 F, its 
destruction by proton capture is rather well known at 
nova temperatures and should not affect the analysis of 
Sect. |. Below «4xl0 8 K, the 17 F(p,7) 18 Ne reaction 
rate is dominated by the effect of direct capture on 
bound 18 Ne levels. This contribution is expected to 
suffer little uncertainty since it has been calculated 
(Wiescher & Kettner 1982; parcia et al. 1991 ) using ex- 
perimental neutron spectroscopic factors from the mirror 
nucleus 18 O. The lowest resonances are located at E r = 



0.586, 0.655 and 0.600 MeV ([Bardayan et al. 1999|) too 
high to contribute, even by their tails, to the rate below 
?s4x 10 8 K. Partial widths are reasonably known through 
experiment s, shell model calculations or from the conju- 
gate levels (Wiescher fc Kettner 1982 Garcia et al. 1991) 



and are such that T p 3> L-y (i.e. ury — luT^) . A 



long standing uncertainty (Wiescher fe Kettner 1982 



Wiescher, Gorres & Thielemann 1988 



Garcia et al. 1991) concerned the location of the £=0, 



0.600 MeV resonance associated with the 4.524 MeV, 3+ 
level. Fortunately, the resonance energy first measured by 
Garcia et al. (1991), has been experimentally confirmed 
very recently by Bardayan et al. (1999). Hence the un- 
certainty associated with the 17 F(p,7) 18 Ne rate appears 
negligible within the domain of nova nucleosynthesis. 

According to NACRE ( [Angulo et al. 1999 1), the 
17 0(p,a) 14 N (Fig. [If) and 17 0(p, 7 ) 18 F (Fig. [$) reac- 
tion rates present large uncertainties at temperatures 
below a few 10 s K. The increase in the rates around 
50xl0 6 (Angulo et al. 1999) is due to the contribu- 
tion of a E r — 66 keV resonance (Landre et al. 1989; 
Blackmon et al. 1995), not included in Caughlan & 
Fowler (1988). The large uncertainty around 2xl0 8 K 
(Angulo et al. 1999), well within the range of tempera- 



ture reached in nova outbursts, is due to the unknown 
contribution of a 179.5 keV resonance associated with 
the = 5786 keV ( 18 F) level. The upper limits for the 
strengths [(p,7) and (p,a)] of this resonance are calculated 
( Angulo et al. 1999| ) from the partial widths extracted by 
Rolfs, Berka & Azuma (1973) and the upper limit for the 
proton width from Landre et al. (1989). As usual, for the 
calculation of the contribution of this resonance to the 
NACRE low, recommended and high rates, its strength 
has been taken equal to 0, 10% and 100% of the ex- 
perimental upper limit. Accordingly, the NACRE recom- 
mended rate is somewhat arbitrary in the temperature 
domain affected by the 179.5 keV resonance (Figs. [I| and 
Figs. ||) 



5. Results: effect of new rates on 18 F synthesis 

New models of a 1.25 M Q ONe nova have been computed 
with the SHIVA code, in order to analyze the effect of 
the new nuclear reaction rates. For the 18 F+p rates, we 
have adopted the low, high and nominal prescriptions de- 
scribed in section 3, whereas for the 17 0+p ones we have 



adopted the NACRE rates (see Table 0). It is impor- 
tant to stress that the energetics of the explosion remains 
practically unchanged with the new 18 F+p and 17 0+p 
reaction rates, since these reactions are not the ones af- 
fecting most the evolution (see Jose & Hernanz (1998) 
and Jose et al. (1999) for a deep analysis of the reac- 
tion rates which have the largest influence in nova ex- 
plosions). Therefore, the maximum temperature attained 
in the burning shell at the base of the accreted enve- 
lope, T ma2 .=2.51 x 10 8 K, the mean kinetic energy and 
the ejected mass, AAf ej - ec =1.8xlO~ 5 M , are unchanged 
with respect to previous models (i.e., Jose et al. (1999)). 

Two previously computed evolutionary sequences, 
with the same white dwarf mass and composition, are 
adopted as old models for the purpose of comparison: a 
case with Wiescher & Kettner (1982) rates (see Jose & 
Hernanz (1998) for the yields and Gomez-Gomar et al. 
(1998) for the 7-ray spectra) and a case with Utku et al. 
(1998) rates (see Hernanz et al. (1999) for the yields and 
the spectra). In both of them, the Caughlan & Fowler 
(1988) updated by Landre et al. (1989) rates for 17 0+p 
were adopted. In Table [l] we show the 18 F yields obtained 
with the new rates, as well as those with the old mod- 
els. We see that yields with the nominal 18 F+p rates are 
smaller by a factor of 30 than those with the older rates, 
as expected. In fact, the reduction in 18 F production is 
even larger (by a factor of 3) than the one obtained with 
Utku et al. (1998) rates (see Hernanz et al. (1999)), as a 
consequence of the inclusion of the tail of the resonance 
at 38 keV (see section 3). 

Concerning the effect of the new 17 0+p rates, a re- 
duction by an extra factor of 2 is obtained when the 
NACRE recommended rates, instead of the old Caughlan 
& Fowler (1988) and Landre et al. (1989) ones, are adopted 
(together with 18 F+p nominal rates; see Table [l]). In 
summary, our nominal yield of 18 F is 4.84xl0~ 5 , which 
is 60 times smaller than the yie ld with the old Wi- 
escher fc Kettner (1982 ) rates ( Jose fc Hernanz 1998 ; 
|G6incz-Gomar et al. 1998 ) and 6 times smaller tha n that 



with Utku et al. (1998) rates ( [Hernanz et al. 1999| ). The 
consequences of the reduced 18 F yields are quite large 
for the early gamma-ray output, which is directly related 
to the amount of 18 F synthesized. The reduced content 
of 18 F in an expanding envelope with similar physical 
and dynamical properties, translates directly into a re- 
duced positron annihilation gamma-ray flux, roughly by 
the same factor as the 18 F decrease. This was shown in 
Hernanz et al. (1999), where complete computations of 
the gamma-ray spectra were performed with the 18 F yields 
obtained with Utku et al. (1998) rates, for different masses 
and compositions of the white dwarf. With the rates pre- 
sented here, the reduction of the predicted fluxes becomes 
larger by a factor of 6, for the combined nominal 18 F+p 
and recommended 1 O+p rates. In summary, if we com- 
pare with results in Gomez-Gomar et al. (1998), where 18 F 
yields were computed with Wiescher & Kettner (1982), the 
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511 keV and continuum fluxes for a 1.25 M Q ONe nova 
should be reduced by a factor of 60. 

Another important point to stress is the impact of 
the uncertainty in the rates on the computed yields. The 
above nominal yields rely upon recommended or nomi- 
nal rates which contain assumed values (^=0.01 or 0.1 
for 18 F+p and 10% of an experimental upper limit for 
17 0+p). Hence, the nominal yield should not be disso- 
ciated from the large remaining uncertainties. From our 
results (see Table [l]) we conclude that the range between 
high and low presciptions for the 18 F+p rates translates 
into an uncertainty (maximum versus minimum yield ra- 
tio) of a factor of 310, a quite large value, whereas for the 
17 0+p rates the uncertainty is a factor of 10. (It is worth 
noticing that the results obtained with the Utku et al. 
(1998) rates fall within the limits obtained with the new 
rates; see Table [l]). The corresponding uncertainties in 
the gamma-ray fluxes are of the same order of magnitude. 
This points out the interest of more accurate determina- 
tions of the 18 F+p and 17 0+p rates. 

6. Conclusions 

We have investigated the 18 F formation and destruction in 
nova outbursts, identified the key reactions (proton cap- 
ture on 18 F and 17 0) and analysed their rates. The proton 
capture rates on 18 F are higher than the Wiescher & Ket- 
tner (1982) ones at nova temperatures due to the contribu- 
tion of the tail of the 659 keV reson a nce whose large mea- 
sured width ( Coszach et al. 1995 ;_ Graulich et al. 1997 ) 



has been indirectly confirmed ( Butt et al. 1998| ). Another 
important contribution comes from the tail of the 38 keV 
resonance which was neglected in previous studies. Its 
nominal contribution is larger than the 659 keV one but is 
proportional to its assumed reduced proton width 9^. The 
strengths (proton widths) of the low lying resonances are 
unknown and induce large uncertainties (factors of 100 to 
1000) in the rates at nova temperatures. We have provided 
updated nominal rates for the two capture reactions to- 
gether with upper and lower limits. The 17 0+p rates also 
display some large uncertainties at nova temperatures ac- 
cording to the recent compilation of Angulo et al. (1999). 

We have used these new nuclear physics results in a 
fully hydrodynamical nova code to calculate the 18 F yields 
in novae for different rates : our low, nominal and high 
18 F+p rates, and the low, recommended and high 17 0+p 
rates from Angulo et al. (1999). These results have been 
compared with models computed with the old Wiescher 
& Kettner (1982) rates ( Jose fc Hernanz 1998[ ), and with 
more recent models with the Utku et al. (1998) rates 
( Hernanz et al. 1999| ). Two important results have been 
obtained. First, there is always a reduction of the amount 
of 18 F synthesized in a nova explosion, with the nominal 
rates for the 18 F+p reactions both alone and combined 
with the recommended rates for the 17 0+p reactions. The 
nominal 18 F yield (nominal 18 F+p and recommended for 



17 0+p) is 4.84 xlO -5 by mass, which is 60 times smaller 
than the one obtained with Wiescher & Kettner (1982) 
rates and 6 times smaller than the one with Utku et al. 
(1998) rates (and Caughlan & Fowler (1988) and Landre et 
al. (1989) for 17 0+p). The impact on the early gamma- 
ray spectrum of the nova is a reduction of the flux by 
the same amount (with respect to Gomez-Gomar et al. 
(1998) and Hernanz et al. (1999), respectively). Second, 
the yields are found to be very sensitive to the rates with 
resulting combined ( 18 F+p and 17 0+p) uncertainties of 
more than three orders of magnitude. This supports the 
need of new experimental and theoretical studies to im- 
prove the knowledge of the 18 F+p and 17 0+p rates and, 
consequently, allow for a larger reliability of the predic- 
tions of annihilation gamma-ray fluxes from novae, to be 
observed by current and future instruments. 
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Fig. .1. Hot CNO cycle. Dashed and dotted circles represent beta unstable and proton unbound nuclei respectively. 
Dashed arrows represent reactions of negligible influence in novae. 

Appendix: Updated rates 

We give here the tabulated reaction rates for proton capture on 18 F resulting from numerical integrations (see § |^ 
and Figs. U and 0). For convenience, we also provide here formulas that approximate the nominal rates (As usual, 
T9 is the temperature in GK, T9XY=T9**(X/Y) and T9LN=LOG(T9) in FORTRAN notations.) 

The following 18 F(p,a) 15 nominal rate formula includes the contribution of the two measured resonances (330 and 
659 keV) with additional contributions from the possible resonances at 38 keV, (with high energy tail) and 287 keV. 

SVDir = 9.13el0/T923*EXP(-18.052/T913-0.672*T92)*(l.+! Tail 
k 0.0231*T913+6. 12*T923+0 . 988*T9+9 . 92*T943+4 . 07*T953) ! 659 keV 
k +5.78e5/T932*EXP(-3.830/T9)+9.91e8/T932*EXP(-7.648/T9)+ ! 330, 659 keV 
k 2.81e-6/T932*EXP(-0.441/T9)+2706.*EXP(3.8319*T9LN-1.3450* ! 38 keV 
k T9LN**2-.0001/T9**3)+4.46e4/T932*EXP(-3.331/T9) ! 287 keV 

The following 18 F(p,7) 19 Ne nominal rate formula contains contributions from direct capture and from the 38, 287, 
330 and 659 keV resonances. 

SVDir = 3.98e7/T923*EXP(-18.052/T913)* ! DC 
k (1 . +0 . 0231*T913+0 . 0885*T923+0 . 0143*T9)+ 
k 1 . 34e3/T932*EXP (-3 . 830/T9) + 1 . 51e4/T932*EXP (-7 . 648/T9) 
k +EXP (-0 . 56781+3 . 6850*T9LN-1 . 3636*T9LN**2- . 0001/T9**3) 
k +8 . 26e-10/T932*EXP (-0 . 441/T9) +12 . 2/T932*EXP (-3 . 331/T9) 
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Fig. .2. Mass fractions of isotopes in the expanding envelope at various stages of the outburst. The horizontal axis 
corresponds to a relative mass coordinate where the origin is set at the surface of the envelope and where the arrow 
shows the limit of the ejected material. The panels are numbered from top to bottom and from left to right and 
correspond to temperatures of approximately 5xl0 7 K, 7xl0 7 K, 10 8 K, 2xl0 8 K, T max =2.5x 10 8 K, and to the last 
phases of the evolution, when the nova envelope has already expanded to a size of R w d ~ 10 9 , 10 10 and 10 12 cm, 
respectively. 
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Fig. .3. Level scheme for 19 F and 19 Ne. Dotted lines relate conjugate levels suggested by Utku et al. 1998. The arrow 
represents the energy range where the presence of resonances (or tails of resonances) affects the reaction rates in the 
domain of temperature corresponding to nova outbursts. The dash-dotted lines represents new experimental data since 
Wiescher & Kettner (1982). 

Table .1. Yields (in mass fraction) of 18 F 1 hour after T max , for a 1.25 M ONe nova, for different prescriptions for 
the 18 F+p and 17 0+p nuclear reaction rates. 



Rates for IS F+p 


Rates for L '0+p 


Reference 


18p 




WK82 


CF88+Lan89 


JH98 


2.80 x 10" 


a 


Utku98 


CF88+Lan89 


HJC99 


2.50 x 10" 


4 


Low 


CF88+Lan89 


This work 


8.19 x 10" 


4 


Nominal 


CF88+Lan89 


This work 


9.78 x 10" 


5 


High 


CF88+Lan89 


This work 


2.98 x 10" 


6 


Nominal 


NACRE low 


This work 


1.58 x 10" 


5 


Nominal 


NACRE rec. 


This work 


4.84 x 10" 


5 


Nominal 


NACRE high 


This work 


1.67 x 10" 


4 



References f or the rates: WK82 ( [Wiescher fc Kettner 1982] ), CF88 (|Caughlari fc Fowler 1988| ), Lan 89 (|Landre et al. 198£| ), 
NACRE ( |Angulo et al. 1999| ). References for the models: JH98 QJose fc Hernanz 1998[ )7TLJC99 ( [Hernanz et al. 1999[ ). 
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E R ™ (MeV) 

Fig. .4. Astrophysical S-factor for the 18 F(p,a) 15 reaction. The solid lines represent contributions from resonances 
whose strengths have been directly measured. Dashed curves represent assumed contributions from resonances corre- 
sponding to known 19 Ne levels. Possible contributions associated to expected levels in 19 Ne are not shown. Note the 
importance of the tail of the i? r =38 keV (E X =6AA9 MeV, 3/2 + ) resonance not considered in Utku et al. (1998). 




E R ™ (MeV) 



Fig. .5. Same as Fig. [^but for the 18 F(p,7) 19 Ne reaction. 
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Fig. .6. 18 F(p,a) 15 rates compared to WK82. The solid lines repr esent the low, no minal and high rates and the 
hatched area the nuclear uncertainty. The dashed line represents the Utku et al. 1998 rate. All rates have been nor- 
malized to the WK82 one. 
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Fig. .7. Same as Fig. [6| but for the F(p,7) 9 Ne reaction 
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Fig. .8. Comparison between NACRE (low, recommended and high) and CF88 rates for the 7 Q(p,a) N reaction 



The origins for the differences are labeled with the energy of the corresponding resonances. (See Angulo et al. 1999 
for more detail.) 




Fig. .9. Same as Fig. [1| but for the 17 0(p,7) 18 F reaction. 
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Table .2. 18 F(p, 7 ) 19 Ne and 18 F(p,a) 15 rates 
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0.623E+06 


1 


250 


0.649E+02 


0.679E+02 


0.982E+02 


0.939E+06 


0.123E+07 


0.189E+07 


1 


500 


0.101E+03 


0.107E+03 


0.169E+03 


0.196E+07 


0.258E+07 


0.394E+07 


1 


750 


0.140E+03 


0.151E+03 


0.259E+03 


0.321E+07 


0.424E+07 


0.655E+07 


2 


000 


0.180E+03 


0.197E+03 


0.361E+03 


0.454E+07 


0.601E+07 


0.948E+07 


2 


500 


0.262E+03 


0.292E+03 


0.581E+03 


0.703E+07 


0.936E+07 


0.155E+08 


3 


000 


0.347E+03 


0.390E+03 


0.800E+03 


0.897E+07 


0.120E+08 


0.209E+08 



